Beyond Concordance Cosmology with Magnification of Gravitational- Wave Standard Sirens 
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We show how future gravitational-wave detectors would be able of discriminating between the 
concordance ACDM cosmological model and up-to-date competing alternatives, e.g. dynamical 
dark energy models (DE) or modified gravity theories (MG). Our method consists in using the 
weak-Iensing magnification effect that affects a standard-siren signal because of its travelling trough 
the Universe's large-scale structure. As a demonstration, we present constraints on DE and MG 
from proposed gravitational- wave detectors, namely ET and DECIGO/BBO. 
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Introduction. — In the last decades, the A cold dark matter 
(ACDM) paradigm has been established as the concordance 
cosmological model. However, the origin of the cosmological 
constant A is still unknown and cosmologists have recently 
proposed several models alternative to it. These models 
mainly rely cither on the introduction of additional scalar or 
vector fields in the Universe's content or on a modification of 
the law of gravity; they thus attempt in finding an agreement 
at least as good as that of ACDM with current cosmological 
datasets, for instance the temperature anisotropy pattern of 
the cosmic microwave background radiation [l[ , the dynam- 
ics of the large-scale structure of the Universe 0, and the 
present-day cosmic accelerated expansion 

Amongst the most-studied alternative models with addi- 
tional fields there is dynamical dark energy (or quintessence; 
hereafter, DE) [3|. Conversely, the class of modified gravity 
(MG) includes a variety of approaches. Actually, the effort 
of modifying the action of gravity dates back to just few 
years after Einstein's seminal papers [e.g. [H, for a historical 
review]. To give some well-studied and representative ex- 
ample, we refer to comprehensive reviews [6|-l8|. It is worth 
noting that a completely different route might be taken, i.e. 
that of trying to mimic (at least a fraction of)A effects with 
the inhomogeneities of the cosmic geometry [9| or backreac- 
tion effect of matter inhomogeneities [l^. This is indeed an 
appealing possibility that needs no new energy component 
of Universe nor modification of gravity. 

In this Letter, we address the challenging and contem- 
porary issue of discriminating amongst all these plausible 
cosmological models by extending the powerful method in- 
troduced in Ref. [Il| . We make use of the weak-lensing mag- 
nification effect on a GW from a compact binary object — 
which is often referred to as "standard siren." It provides a 
unique way to measure the luminosity distance to the source 
(l^ . Proposed ground- and space-based GW detectors, such 
as Einstein Telescope (ET) [13[, DECThertz interferometer 
Gravitational- wave Observatory (DECIGO) [3], and Big- 
Bang Observer (BBO) [lH, will detect a large number of 
neutron-star (NS) binaries, hence allowing us toprobe the 
cosmic expansion with unprecedented precision [ll, 16- 1^. 

Ref. [m showed that the weak-lensing magnification ef- 



fect on a GW, instead of being a noise term, in fact carries 
useful information about density fluctuations. Specifically, 
they demonstrated that future GW detectors will be able to 
measure a lensing convergence power spectrum. It has also 
been shown that this method can be used to constrain cos- 
mological parameters by taking the cosmological-paramctcr 
dependences of the lensing magnification distribution into 
account 2^ 21 1. Here, we aim to step forward their method 
by treating the angular power spectrum of the lensing mag- 
nification as a signal and use the method for cosmological 
parameter estimation and model selection between DE and 
MG [e.g. |22|. It is in a sense a higher- level question than pa- 
rameter estimation. In parameter estimation, one assumes 
a theoretical model within which one interprets the data, 
whereas in model selection, one wants to know which theo- 
retical framework is preferred, given the data (regardless of 
the parameter values). 

GW standard siren. — For a single binary system, the 
Fourier transform of the GW waveform is expressed as a 
function of frequency v [23| , 
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where c^l is the luminosity distance, = (1 + z)Mc is 
the redshifted chirp mass, and Mc is the proper chirp mass. 
The constant A is given by ^ = l/(-\/67r^/'^), which includes 
the factor ^4/5 for a geometrical average over the inclina- 
tion angle of a binary. The function '^{v) represents the 
frequency-dependent phase including terms up to the or- 
der of the restricted 1.5 post-Newtonian (PN) approximation 
[2^. For GW sources at cosmological distance, the phase is 
modulated due to cosmological acceleration. However, its 
contribution is small and can be ignored if we assume that 
source redshifts are obtained by electromagnetic (EM) ob- 
servations [igj . 

GW observations measure the mass parameter = 
(1 + z)Mc and cannot provide source redshift separately. 
Therefore, the redshift has to be determined via EM ob- 
servations. There are a few proposals for standard sirens 
which can lift the redshift degeneracy 13, 13, 25 1. However, 
to clarify the feasibility of these methods, further studies 
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FIG. 1. Accuracy on the luminosity distance measurement with 
DECIGO (red) and ET (green) for 3-year observation. 



are necessary, since each method is based on different prior 
knowledge on the GW sources or different detector sensitiv- 
ities. On the contrary, in this Letter we assume that source 
redshifts are known by the EM follow-up observations of host 
galaxies or the EM counterparts of binary merger events. 
Since the assumption that redshifts of all sources are identi- 
fied is too optimistic to be justified, we consider the fraction 
of the redshift identification amongst all GW sources to be 
10-3, 10-2, 10-1 and 1. 

The distance accuracy AdL{z)/dL{z) for DECIGO in a 
ACDM Universe has been recently estimated via Fisher ma- 
trix analysis [l^. Here, we follow the same approach to also 
obtain what will be expected for ET, using the noise curve 
given in [l^. The results that will be achieved for 3-yr ob- 
servation time of a single source is depicted in Fig. [1] 

The number of NS binaries in the redshift interval [z, z + 
dz] observed during the observation time Tobs is given by 
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where x{^) = dL{z)/{\ + z) is the radial comoving distance 
and h{z) = nos{z) with s{z) 1 + 2z for z < 1, |(5 — z) 
for 1 < z < 5, and for z > 5. Although the normalisation 
of n is still uncertain, we adopt the most recent estimate 

no = 10-VMpcVyr 

Cosmology beyond the concordance model. — From struc- 
ture formation, we know that the A-like component cannot 
clump on the scales covered by the galaxy surveys and below. 
Thus, it must be something different from standard matter 
and radiation and even from dark matter. To fi nd p ossible 
solutions to the "cosmological constant problem" [29| , many 
routes have been tested. Probably, the most well-known is 
scalar- field DE It is known that a self- interacting scalar 



field with negligible kinetic energy behaves like a cosmologi- 
cal constant. The main question about A is why it is not zero 
if it is so incredibly small. DE provides a fairly satisfying an- 
swer to this with the idea of a time-dependent cosmological 
constant, whose present-day value matches our observations. 
Here, we use the Chevallier-Polarski-Linder parameterisa- 
tion 3^, 31 1 , where the equation-of-state parameter of DE 
is a Taylor expansion around a cosmological constant with 
coefficients wq and Wa, viz. wue{z) ~ wo + z/{z + l)wa- The 
ACDM limit is reached when wq —1 and Wa — > 0, these 
arc therefore the values we set as fiducial. 

Regarding MG, we follow a phenomenological approach. 
We assume that MG can mimic the ACDM expansion his- 
tory. Nonetheless, it is known that MG alters the growth of 
cosmic structures. According to the literature, we therefore 
introduce two free functions ^(fc, z) and r]{k, z) which nat- 
urally appear in the perturbed field equations and encode, 
respectively, any modification to the Poisson equation and 
to the ratio of the metric potentials $ and 5*, i.e. 
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where A is the gauge-invariant comoving density contrast 
of background energy density p. General relativity implies 
IJ,{k,z) = ij{k,z) = 1. whereas in MG their functional form 
depends on the model adopted. Since we are here interested 
to study the capabilities of gravitational-wave detectors in 
detecting any MG signature — rather than scrutinise specific 
MG theories — we follow Ref. [32j, who assume the approxi- 
mate expressions 
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thence, = ^J■{z ~ 0) and 770 = ?/(z = 0) are the MG 
parameters. Since we are interested in testing gravity at late 
times, we have that /i(z) = r]{z) = 1 for z ^ Zs, where Zs is 
a threshold redshift and Az is the width of the transition. 
As fiducial values, we take the ACDM limit, viz. {/xq, 770} = 
{1, 1}. The other parameters are kept fixed at z^ = 2 and 
Az = 0.05 [33. 

It is worth giving a final remark on the evolution of mat- 
ter fluctuations. Unlike the linear growth of perturbations, 
which can be described analytically, the non-linear regime 
has to be explored numerically. At the time being, there is 
not a thorough agreement in the literature on how DE and 
MG behave in the non-linear regime. However, all the at- 
tempts hitherto made point towards a modification of the 
HALOFIT procedure [33[. We therefore proceed as follows: 
we adopt halofit as non-linear procedure, but shall not 
consider angular wavenumbers i > 1000, to stay safely out 
of the highly non-linear regime. 

Method. — Following the approach outlined in Ref. [ll[ , we 
use the wcak-lcnsing magnification power spectrum as a car- 
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dinal observable. The underlying idea is straightforward: 
in the geometrical optics approximation, GWs follow nuU- 
geodesics of the spacetime. Hence, their propagation from 
the source towards us is formally the same as that of light 
rays. As a consequence, such signal will suffer distortions due 
to the distribution of matter in the Universe in its passing 
trough it. In other words, the potential wells of the cosmic 
large-scale structure will bend the GW trajectories in the 
same way they curve photon paths. For the well-developed 
theory of gravitational lensing, we refer to the literature [3J] . 

Here, we choose to exploit the power of tomography, where 
the redshift distribution of the sources — NS in the present 
case — is further subdivided into redshift bins. This is a sig- 
nificant improvement regarding the present technique ap- 
plied to GWs. Indeed, this will allow us to elucidate the time 
dependence of DE and MG. The magnification tomographic 
matrix C^^ {£) as a function of the angular wave-number £ 
reads 
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is the weak-lensing selection function in the i-th redshift 
bin; zh the survey depth; dN^^g/dziz) the redshift distri- 
bution of NS in that bin, with N^g its normalisation, viz. 

JdzdN^g/dziz); Hq the Hubble constant; the total mat- 
ter fraction; and P^ik, z) the power spectrum of matter fluc- 
tuations as a function of the physical scale k and the redshift 
z. Here, wc use Limber's approximation, which sets k ~ £/x 

We perform a Fisher matrix analysis [36| for the set of 
cosmological parameters ??q, = {f2m, n^, erg, wq, Wa} for DE 
and {flm, Us, as, Mo, '7o} for MG, where is the slope of 
the primordial power spectrum and ag is the rms mass fluc- 
tuations on a scale of 8h~^ Mpc; the extra-ACDM parame- 
ters are {wq, Wa} and {/io, rjo} for DE and MG, respectively. 
The Fisher matrix elements can be written as 
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where /sky is the fraction of the sky covered by the survey 
under analysis, 
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NS 



is the covariance matrix of the signal, Sij the Kronecker sym- 
bol and a^{zi) the instrumental error on the magnification 
estimate as a function of the redshift bin. We will spend 
more comments on it in the following. 



As we explained, a GW detector itself will not mea- 
sure source redshifts. Therefore, we have to rely on cross- 
identifications made by ancillary surveys. To this purpose, 
there are many forthcoming spectroscopic and photomet- 
ric experiments that will provide large-scale catalogues of 
galaxies. To enumerate some of them: Pan-STARRS 
SkyMapper S, SDSS-III d!], Euchcfl H^, the LSST 
or DES [42]. Such a large number of surveys covering the 
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whole sky ensures that a non-negligible fraction of the NS de- 
tected by ET/DECIGO may in principle be cross-identified. 
This enables us to perform the tomography. Nevertheless, 
though one of the main features of the dAf]^*g/dz(z) is its 
depth, none of the surveys considered will reach objects at 
z > 2. Hence, we fix the survey depth of Eq. ([U to z// = 2, 
and we further subdivide the redshift distribution of sources 
into five equally-spaced bins. 

Regarding the noise term ct^^ in Eq. pO)) . it is related to 
the instrumental error on the estimate of the luminosity dis- 
tance dLiz) from the GW signal in Fig.[T] The instrumental 
scatter is then simply given by cr^ = H^di^j d^- Since we are 
calculating the magnification tomographic matrix, we need 
to define 17^^(2:^), namely as a function of the redshift bin. 
Here, we proceed by averaging over the extension of the bin. 
Moreover, we weight the signal with the redshift distribution 
of NS. Whence, 
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Discussion of the results. — To show the main results of 
our analysis, we make use of the Figure of Merit (FoM) for 
the DE and MG parameters. It has been introduced by 



the Dark Energy Task Force [43|, and it is proportional to 
the area encompassed by the ellipse representing the 68.5% 
confidence level in the two-parameter plane. Specifically, for 
the parameter pair 'da-'&p, it reads 
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det (F- 
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In Fig.[2]we present FoM(wo, Wa) and FoM(/^o, %), the FoMs 
for the DE and MG extra-ACDM parameters, versus the 
fraction of ET/DECIGO sources cross-identified by other 
surveys. As a comparison, we also show our results in com- 
binations with priors from the Planck satellite |44|. 

As expected, the FoMs present a direct proportionality 
with the fraction of detected sources. Similarly, constraints 
from DECIGO are always more stringent than what obtained 
from ET — the accuracy on the luminosity distance measure- 
ments of the latter experiment is worse than that of the 
former. CMB priors from Planck are definitely dominant 
when the number of detected sources is small, whilst for in- 
termediate fractions they help in enhancing the FoMs from 
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FIG. 2. FoM(u'o, Wa) (left panel) and FoM(/^o, 570) (right panel) as a function of the fraction of detected sources for DECIGO (squares, 
red) and ET (diamonds, green). Solid lines refer to GW detectors only, whilst dashed lines show the results in combination with 
Planck priors. 



TABLE I. Forecasted la errors on the extra-ACDM parameters 
for 100% source-redshift identification. 





wo 


Wa 


MO 


'70 


DECIGO 


0.022 


0.062 


0.032 


0.23 


ET 


0.098 


0.42 


0.15 


0.78 



DECIGO and ET. If all the sources were cross-identified, 
the larger contribution would be that of GW detectors. We 
summarise our more promising results in Table HI 

Besides, it is easy to appreciate how powerful our method 
is in constraining DE. This is mainly because DE alters the 
background evolution of the Universe, leading to the present- 
day accelerated expansion. Therefore, the magnification 
signal due to GWs receives a contribution from both the 
background level and the perturbations to the matter power 
spectrum. On the other hand, since in the phenomcnologi- 
cal model for MG we use here the Hubble expansion H{z) 
is exactly the same as in the ACDM model, FoM(/xo, 770) 
depicts only the sensitivity of GW magnification from the 
regime of perturbations. However, it is worth reminding the 
reader that a specific MG model would also modify the Uni- 
verse's background evolution; this means that the method 
we present here is even more effective, when applied to a 
specific theoretical model. 

As another possibility, there are inhomogeneous models 
that mimic ACDM 0, [l^l ■ In principle, our method may 
be used in this case too because GW observation is sensi- 
tive to both background and perturbative evolution. The 
application to rcdshift drift in the Lcmaitre-Tolman-Bondi 
void model has been investigated in Ref. [Zij. However, the 



description of cosmological perturbations is trickiesome in 
these models, and goes beyond the purpose of the present 
work. 

Finally, we note that the method in Refs [2C 



21| is based 



on ID (angular-averaged) distribution of lensing magnifi- 
cation. On the other hand, our method relies on angu- 
lar correlation of magnification and is qualitatively different 
from theirs. Thus, the comparison of the sensitivity is not 
straightforward at this stage, but it would be interesting to 
compare these two methods in the future work. 

Conclusions. — In this Letter, we have shown that the fu- 
ture gravitational- wave detectors, ET and DECIGO/BBO, 
will be able to probe for matter inhomogeneities of the Uni- 
verse through lensing magnification and bring us crucial in- 
formation beyond background-level cosmology. We want to 
emphasise that, although this idea was already proposed fll| . 
we exploit it here to address a more fundamental question, 
namely if our current concordance cosmological model were 
to be preferred by future data, or competing models such 
as DE and MG represented an actual, more viable alterna- 
tive. From the starting point presented in the present work, 
further model-selection anal yse s may be implemented, e.g. 
the use of the Bayes factor [22] ■ However, we find it more 
straightforward and immediate to show our results to the 
two communities of GWs and cosmology by exploiting the 
widely-known FoMs. 

Finally, we note that GW observations are completely in- 
dependent of EM observations and allow us to cross-check 
the methods used in the EM observations, e.g. distance cal- 
ibration of type la supernova. For the feasibility of GWs 
as a cosmological probe, it is important to enhance the frac- 
tion of source rcdshift identifications by further investigating 
the EM counterparts of NS binary merger and developing a 



5 



large-scale follow-up campaign dedicated for GW events. 
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